Exchange interactions and Curie temperature in (GaMn)As 
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We use supercell and frozen-magnon approaches to study the dependence of the magnetic interac- 
tions in (Ga,Mn)As on the Mn concentration. We report the parameters of the exchange interaction 
between Mn spins and the estimates of the Curie temperature within the mean-field and random- 
phase approximations. In agreement with experiment we obtain a nonmonotonous dependence of 
the Curie temperature on the Mn concentration. We estimate the dependence of the Curie temper- 
ature on the concentration of the carries in the system and show that the decrease of the number 
of holes in the valence band leads to fast decrease of the Curie temperature. We show that the 
hole states of the valence band are more efficient in mediating the exchange interaction between Mn 
spins than the electron states of the conduction band. 



I. INTRODUCTION 

An important current problem on the way to the prac- 
tical use of the spin-transport in semiconductor devices is 
the design of the materials that make possible the injec- 
tion of the spin-polarized electrons into semiconductor at 
room temperature. One of the promising classes of ma- 
terials are the diluted magnetic semiconductors (DMS) 
of the III-V type. A strong interest to these systems 
were attracted by the observation of the ferromagnetism 
in (GaMn)In0 and (GaMn)Afi with the Curie temper- 
ature of Gao.947Mno.o53As as high as 110 K. To design 
the material with Curie temperature higher than room 
temperature the knowledge of physical mechanisms gov- 
erning the exchange interactions in these systems is of 
primary importance. 

The theoretical works on the ferromagnetism in the 
DMS systems can be separated into two groups. First 
group models the problem with an effective Hamiltonian 
containing experimentally determined paramctcina. This 
part of the studies is recently reviewed in Rcfs.0cl. This 
paper belongs to the second group of studies that are 
based on the parameter-free calculations within the den- 
sity functional theory. Several calculations iw-ffi, been 
performed recently along this line. In RefsJjoEl, the 
coherent potential approximation (CPA) was used to 
study the magnetic structure, density of states, toteln 
ergy, and chemical trends in (III,Mn)V. In Refs.t 
the calculations were performed for a series of magnetic 
states of (III,Mn)V using supercells of the zinc-blende 
structure and focusing basically on the same physical 
quantities as the CPA studies. Since the CPA is a single- 
site theory which neglects any short range order in the 
disordered subsystems the Qir/Wf'^'^ supercell approaches 
are complimentary. In Refs.l3l3'E2l, the influence of the an- 
tisite defects on the magnetic properties of the DMS is 
discussed. In Ref.llil, the total energy of various coUinear 
magnetic configurations are used to estimate the param- 
eters of the exchange interaction between the 3d atoms 
forming clusters in III-V semiconductors. 

The purpose of the present paper is parameter-free cal- 
culation of the exchange interactions and Curie tempera- 
ture in (GaMn) As for various concentrations of Mn. The 



study is based on the supercell approach. 



Density functional theory (DFT) has proved to be 
very successful in the parameter-free description of the 
ground state magnetic properties of complex systems 
(see, e.g.E3). Recently much attention has been devoted 
to the application of the methods of the DFT to the 
studies of the low-energy excitatior 
and the magnetic phase transitional^ 
The Stoner theory, which relates the finite temperature 
effects to the temperature variation of the Fermi-Dirac 
distribution appeared to be unable to describe the tem- 
perature propepties of the magnetic systems with itin- 
erant electronsBj A necessary feature of the theoretical 
description of the finite temperature effects is an account 
for transversal fluctuations of the local magnetization. A 
most consequent method of the calculation of the low- 
energy magnetic excitations is based on the evaluation of 
the nonuniform antiiequency dependent enhanced mag- 
netic susceptibility.lla This approach is, however, compu- 
tationally very complicated and up to now has been suc- 
cessfully applied to the simplest magnetic systems only. 



A more tractable approach to the study of both the 
spin-wave excitations and the thermodynamics of mag- 
netic systems is based on an. ndiabaJif^paimcnt of the 
atomic magnetic moments.E^EJtjIlj'OoEErEa In this ap- 
proach the account for noncoUinear conflgurations of 
atomic moments is essential. An effective method for 
the estimation of the parameters of the interatomic ex- 
change interaction and spin-wave-energies is suggested 
by the frozen-magnon approach.EsEil This approach is 
based on the total energy calculation for spiral magnetic 
configurations. Because p£ the generalized translational 
periodicity of spin spiral a^n such calculations can be per- 
formed very efficieiitlyO Additional help is provided by 
the force theoremll3 that allows to use band energy of 
non-selfconsistent frozen magnon states for the estima- 
tion of the total-energy differences. In the present paper 
this approach is used to study the exchange interactions 
and Curie temperature in Gai_xMnxAs with various con- 
centrations of the Mn impurities. 
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II. CALCULATIONAL SCHEME 

The calculations are based on the supercell approach 
where one of the Ga atoms in a supercell of zinc-blende 
GaAs is replaced by the Mn atom. The concentration 
of Mn depends on the size of the super cell. The fol- 
lowing concentrations x have been studied in this work: 
0.25, 0.125, 0.0625, and 0.03125. For comparison, the 
calculation of the density of states (DOS) of pure MnAs 
and GaAs are also performed. In oder to better identify 
trends we investigate here a concentration range that is 
much wider the one accessible to experiment sofar. 

The calculations were carriedpOut with the augmented 
spherical waves (ASW) methodBj In all calculations the 
lattice parameter was chosen to be equal to the experi- 
mental lattice parameter of GaAs. Two empty spheres 
per formula unit have been used in the calculations. The 
positions of empty spheres are (0.5, 0.5, 0.5) and (0.75, 
0.75, 0.75). Radii of all atomic spheres were chosen to 
be equal. Depending on the concentration of Mn, the 
super cell is cubic (x=25%, a x a x a, and x=3.125%, 
2a X 2a X 2a) or tetragonal (x=12.5%, a x a x 2a and 
6.25%, 2a X 2a X a). 

To describe the exchange interactions in the system we 
use an effective Heisenberg Hamiltonian of classical spins 



H, 



(1) 



where Jij is an exchange interaction between two Mn 
sites and is the unit vector pointing in the direc- 
tion of the magnetic moments at site i. 

To estimate the parameters of the Mn-Mn exchange 
interaction we performed calculation for the following 
frozen-magnon configurations: 



— const. 



q 



(2) 



where 9i and (pi are the polar and azimuthal angles of 
vector Bi, Ri is the position of the ith Mn atom. The 
directions of the induced moments in the atomic spheres 
of Ga and As and in the empty spheres were kept to be 
parallel to the z axis. 

It can be shown that within the Heisenberg model (|l|) 
the energy of such configurations can be represented in 
the form 

E{e,q)^Eoie)-^-^Jiq) (3) 

where Eq does not depend on q and J(q) is the Fourier 
transform of the parameters of the exchange interaction 
between pairs of Mn atoms: 

>^(q) = E '^OJ exp(iq • Roj). (4) 

In Eq. (H) angle 9 is assumed to be small. Using J(q) one 
can estimate the energies of the spin-wave excitations: 



c.(q) = ^[^(0)-J(q)]^^ 



{Ei6,q)~E{0,O)) 



(5) 



where M is the atomic moment of the Mn atom. Per- 
forming back Fourier transformation we obtain the pa- 
rameters of the exchange interaction between Mn atoms: 



^"J " ^ E exp(-«q • Roj) J(q). 



(6) 



The calculation of E{9, q) for different Mn concentra- 
tions has been performed for uniform meshes in the first 
BZ for 9 = 30°. The symmetry of the system was em- 
ployed to reduce the amount of calculations. For cubic 
supercells (Mn concentrations 25% and 3.125%) the crys- 
tal structure possesses 24 point symmetry operations, for 
the tetragonal unit cells (Mn concentrations 12.5% and 
6.25%) the number of point operations is 8. Since the 
energy as a function of q is invariant with respect to the 
reversal of q the irreducible part of the BZ is |^th and 
jgth correspondingly for cubic and tetragonal supercells. 

The number and type of the exchange parameters de- 
termined in the back Fourier transformation with Eq. (|^) 
are uniquely determined by the q mesh. The parameters 
obtained in this procedure guarantee that J(q) calculated 
according to Eq (^ reproduce exactly the calculated val- 
ues of the total energy. Calculation for a more dense q 
mesh gives additionally the parameters of the exchange 
interaction between more distant atoms. 

The Curie temperature was estimated in the mean-field 
approximation (MFA) 



J#0 



(7) 



and random phase approximation (RPA) 

M 2^ 
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To evaluate the RPA value of the Curie temperature the 
uj{q) was considered continuous in the q space. In a small 
sphere with the center at q = the singular function -^^^ 

was approximated by the function and then replaced 
by a continuous function which has the same value and 
slope at the sphere boundary. The difference between 
the singular and continuous functions was integrated an- 
alytically. The regular function obtained was integrated 
numerically. 

Since the T^^^ and T^^^ are given by the arithmetic 
and harmonic averages of the-spin-wave energies, T^^^ 



is always larger than T^^^£^ Physically this difference 
can be explained by an increased role of the low-energy 
excitations in the case of RPA. In the case of the fer- 
romagnetic 3d metals the RiEA gives, in general, better 
agreement with experiment .E3 (For a recent application 
of the RPA within a model-Hamiltonian approach see, 
e.g., Ref.F 
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FIG. 1: The DOS of Gai_^Mn^As. The DOS is given per 
unit ceU of the zinc-blende crystal structure. The DOS 
above(below) the abscissas axis corresponds to the spin- 
up(down) states. 
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III. DENSITY OF STATES 

We begin with the discussion of the trends in the vari- 
ation of the electron DOS (Fig. |^). For all Mn concen- 
trations studied there is an energy gap in the spin-down 
DOS in an energy interval containing the Fermi level or 
close to it. This gap is about 0.1 Ry in MnAs and de- 
creases to the value of about 0.05 Ry for lower Mn con- 
tent. In MnAs and Gao.75Mno.25As, the Fermi level is 
slightly above the gap. For MnAs, the energy distance 
between the Fermi level and upper edge of the gap is 
about 0.015 Ry. For x — 0.25, this distance is less than 
0.005 Ry. For lower Mn concentrations the Fermi level 
lies within the gap moving from the upper to the lower 
part of the gap with decreasing x. This means that for 
concentrations of Mn less than 25% the calculated ground 
state is half-metallic. This property is very impa£|tant for 
the efficient spin-injection into semiconductor .E3 In the 
half-metallic ferromagnetic state the value of the mag- 
netic moment per supercell is integer (Table |). 

The situation with spin-up DOS is different. In MnAs 
there is no energy gap close to the Fermi level. For 
a;=0.25 a small gap appears at the energy 0.07 Ry above 



TABLE I: Magnetic moments in Gai-^^Mna: As. There are 
shown the Mn moment, the induced moment on the nearest 
As atoms, and the magnetic moment of the unit cell (supercell 
in the case of x 7^ 1). All moments are in units of /is. 

X 

1 0.25 0.125 0.0625 0.03125 
Mn 3.76 3.85 3.88 3.94 3.95 
As -0.18 -0.046 -0.046 -0.036 -0.032 
cell 3.65 3.98 4.00 4.00 4.00 

the Fermi level. With decreasing x, this gap increases 
and its lower edge becomes closer to the Fermi level. At 
the Mn concentration of 6.25% there is an overlap of the 
spin-up and spin-down gaps and an energy gap appears in 
the total DOS. The gap in the total DOS can be treated 
as a gap between the valence and conduction bands of 
the system. The upper part of the valence band is not 
occupied and contains holes. The hole states are of the 
spin-up type. 

Comparison of our DOS with the corresponding DOS 
available in the scientific literature shows good agree- 
ment. Thus, our DOS for the Mn concentratiop of 3.125% 
is close to the corresponding DOS from RefB calculated 
with the use of pseudopotential approach. Also the par- 
tial DOS for x=6.25% presented in Reffl is in good agree- 
ment with the curves obtained in our calculations. 

The values of the calculated moments in the Mn atomic 
sphere are collected in Tableland are close to 4/1^. Simi- 
lar values have been obtained in other calculatiopa within 
the density functional theory (see, e.g., Refs.BEltiJ). Tak- 
ing as an example the system with x = 3.125% we find 
that the contribution of the 3d electrons into the Mn 
moment is 3.83/iB with the rest 0.12/iB coming from the 
4s and 4p electrons. Note that the total number of the 
3d electrons in the Mn sphere is 5.30. The difference 
between the number of the 3d electrons and their contri- 
bution into the spin moment results from the presence of 
0.73 spin-down 3d electrons in the Mn atomic sphere (see 
the spin-down Mn-DOS in Fig. ||). The hybridization be- 
tween the Mn 3d states and the states of the valence band 
of GaAs is crucial for the appearance of the occupied 
spin-down 3d states. In some model-Hamiltonian studies 
a physical picture is used which considers the Mn 3d elec- 
trons as strongly localized and forming an atomic spin of 
5 = |. The density-functional-theory calculations show 
that this picture, although useful in qualitative studies, 
does not take into account some important features of 
the Mn 3d states. 

One of the important issues in the magnetism of the 
DMS is the spatial localization of the hole states. Com- 
parison of the DOS of pure GaAs and the GaMnAs with 
X — 3.125% helps to get insight into the physical mecha- 
nism of the formation of the hole states. The replacement 
of one Ga atom in the supercell of GaAs by a Mn atom 
does not change the number of the spin-down states in 
the valence band. In the spin-up chanel there are, how- 
ever, five additional energy bands which are related to the 
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FIG. 2: The partial Mn-DOS for Gai_ 
given per Mn atom. 



,Mn^As. The DOS is 



FIG. 3: Hole distribution and induced moments (in fis) for 
Gai-xMn^iAs with x — 0.03125%. All values are given for 
the coordination spheres. The numbers of atoms in the coor- 
dination spheres are given at the top of the picture. 




_i , I , lUI 

-0.4 -0.2 



-0.4 -0.2 



E(Ry) 



Mn 3d states. Since there are five extra energy bands and 
only four extra electrons (the atomic configurations of Ga 
and Mn are 4s^4p^ and 3d'^4s^) the valence band is not 
filled and there appear unoccupied (hole) states at the 
top of the valence band. The integrated number of the 
hole states is exactly one hole per Mn atom (correspond- 
ingly, one hole per supercell).The distribution of the hole 
in the supercell for the Mn concentration of 3.125% is 
shown in Fig. ||. About 18.5% of the hole is in the Mn 
sphere, 22.5% in the first coordination sphere of the As 
atoms. Correspondingly about 60% of the hole is outside 
of the first coordination sphere of the impurity. About 
58.3% of the hole is on the As atoms, 17.7% on Ga sites, 
5.5% in the empty spheres. Thus, the hole states are 
rather delocalized. Even the most distant As atom con- 
tains 4.5% of the hole. The delocalized character of the 
hole states is an important factor in mediating the ex- 
change interaction between Mn atoms. 

Another important quantity characterizing the local- 
ization of the valence-band states about impurity is the 
values of the induced moments on various atoms. The 
values of the atomic moments for x — 3.125% are shown 
in Fig.^. It is seen that even at the As atom most distant 
from Mn impurities there is substantial spin polariza- 
tion that provides an efficient exchange path between Mn 
atoms. The dependence of the induced moment on the 
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distance from Mn is nonmonotonous for both As and Ga. 
This behavior can be related to the spin density oscil- 
lation in the Ruderman-Kittel-Kasuya-Yoshida (RKKY) 
theory. The detailed analysis of the formation of the in- 
duced atomic moments shows, however, that the physics 
here is more complex than the physics considered by the 
RKKY theory since an essential role in the formation of 
the induced magnetic moments plays the hybridization of 
the As and Mn states. In particular, the negative sign of 
the moment of the As atoms is explained by the property 
that the empty (hole) states have large As contribution. 
Since these states are of the spin-up type the spin-down 
As states become more occupied than the spin-up ones 
leading to a negative induced moment. 



IV. EXCHANGE PARAMETERS AND CURIE 
TEMPERATURE 

The calculated exchange parameters are presented in 
Fig. g There are a number of conclusions that fol- 
low from the analysis of this figure. First, the Heisen- 
berg model with the interaction between the first near- 
est neighbors only is not able to describe the magnetism 
of the system. For two Mn concentrations {x = 25% 
and X = 3.125%) the first nearest-neighbor interaction is 
even antiferromagnetic, for x ~ 12.5% it is positive but 
small. Second, the exchange interactions is rather quickly 
decreasing with increasing distance between atoms. In 
Fig. we show the variation of the Curie temper- 
ature (M) with increasing number of the contributing 



FIG. 4: The parameters of the exchange interaction between 
Mn atoms (upper panel) and the variation of Tq^^ with 
increasing number of the contributing coordination spheres 
(lower panel). The abscissa gives the radius of the coordina- 
tion sphere in the units of the lattice parameter of the zinc- 
blende crystal structure. 



FIG. 5: Curie temperature of Gai_a;Mn2,As. 
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coordination spheres. For instance, for concentrations 
X — 25%, 12.5%, 6.25% no noticeable contribution to the 
Curie temperature is obtained from the interactions be- 
tween Mn atoms at the distances larger than 3a. Third, 
the dependence of the exchange parameters on the dis- 
tance between Mn atoms is not monotonous. At the dis- 
tance of 2a the interatomic interaction is negative for all 
concentrations considered. Although the theory of the 
RKKY interaction is not sufhcient to describe the mag- 
netism of the system, it is instructive to compare the 
characteristic length of the variation of the calculated 
exchange parameters with the period of the oscillation of 
the exchange parameters which follows from the RKKY 
theory. Since the number of holes is exactly one per su- 
percell the characteristic volume of the Fermi sphere in 
a simple single-band free-electron model is exactly the 
volume of the Brillouin Zone corresponding to the given 
supercell. Thus for x — 25% one can expect the oscil- 
lations with period close to a and for x = 3.125% with 
period close to 2a. Indeed, the analysis of the calculated 
exchange parameters shows that the first and second os- 
cillation of the exchange parameter for x = 25% take 
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place at the distance close to a (from la to 2a and from 
2a to 3a) . On the other hand the characteristic range of 
the variation of the exchange parameter for x = 3.125% 
is larger (from 2a to 4a). Thus there is a correlation be- 
tween the number of holes and the characteristic range 
of the variation of the exchange parameters. 

The calculated Curie temperatures are presented in 
Fig. H Both MFA and RPA approximations give sim- 
ilar dependence of the Curie temperature on the Mn 
concentration: it increases at small x has maximum at 
X = 12.5% and decreases for higher values of x. As 
already stated in Sect. II, the T^^^ is always larger 
than T^j^^. However, even T^^^ exceeds substantially 
the experimental values of the transition temperature. A 
possible explanation of this difference is the presence of 
the As antisites and other donor defects that decrease the 
concentration of the holes-.m-the materials studied exper- 
imentally (see, e.g., Refs.cHZlliy). We address this issue in 
Sect. 0. 

The nonmonotonous dependence of the transition tem- 
perature on the Mn coMjipitration is in qualitative agree- 
ment with experiment .Eal23 This nonmonotonous behav- 
ior is a consequence of the competition between different 
trends arising from the increase of the Mn concentra- 
tion. On the one hand, the amplitude of the effective 
exchange interaction between the Mn moments through 
the valence-band states of GaAs increses with deacreas- 
ing distance between Mn atoms. Also the direct overlap 
of the Mn states increases. These features produce a 
trend to the increase of the Curie temperature with in- 
creasing X. On the other hand, an increased interaction 
of the states of different Mn atoms results in the broad- 
ening of the features of the partial Mn DOS at the Fermi 
level: As is clearly seen in Fig. 0, at a; = 3.125% there is 
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a narrow peak at the Fermi energy that is replaced by a 
broader structure with increasing x. Simultaneously the 
Mn magnetic moment decreases (Table ^ . These prop- 
erties produce the trend to decreasing Curie temperature 
with increasing Mn concentration. 

Interestingly, the induced moment on the nearest As 
atoms is the same for 25% and 12.5% and decreases for 
6.25% and 3.125% although the inducing Mn moment is 
monotonously increasing (Table |). This shows that the 
nearest environment of each Mn atom is influenced by 
other Mn atoms even for the lowest Mn concentration. 



V. DEPENDENCE OF THE CURIE 
TEMPERATURE ON THE NUMBER OF HOLES 

The experimental studies show that the concentration 
of holes in (GcuMn)As is lower than the concentration of 
the Mn atoms Ej One of the important factors leading to 
the low concentration of holes is the presence of the As 
antisites (Asca)- Since As has two more valence electrons 
compared with Ga, each As antisite compensate the holes 
produced by two Mn atoms. 

Here we use a simple rigid-band model to study the de- 
pendence of the Curie temperature on the concentration 
of antisites and other nonmagnetic donors. We assume 
that the electron structure calculated for (Ga,Mn)As 
with a given Mn concentration will basically be preserved 
in the presence of defects. The main difference caused by 
the antisites concerns the occupation of the bands and, 
respectively, the position of the Fermi level. In Fig. ^ 
we show the dependence of the Curie temperature calcu- 
lated in the MFA as a function of the carrier number for 
X — 6.25%. The calculated curve reveals strong depen- 
dence of the Curie temperature on the number of holes 
in the valence band. In agreement with the commonly 
accepted picture of the hole-mediated ferromagnetism in 
(GaMn)As the decrease of the number of holes leads to 
the decrease of the Curie temperature. At the number 
of additional electrons close to 0.85 (the relation of the 
hole concentration with respect to the Mn concentration 
^ = 0.15) the ferromagnetic and antiferromagnetic inter- 
actions compensate and the value of the mean field acting 
on the Mn spins in the ferromagnetic state becomes zero. 
Further decrease of the number of holes leads to negative 
values of the mean field. The negative value of the Curie 
temperature in Fig.^ means that in the assumed ferro- 
magnetic configuration the antiferromagnetic exchange 
interactions prevail resulting in a negative exchange field 
acting on each Mn spin. The ferromagnetic state is un- 
stable. 

At n = 1 the valence band is full. For larger n the 
conduction band is occupied. Discontinuity in the char- 
acter of the occupied states at n = 1 results in a kink in 
the n-dependence of the value of the mean field. With 
occupation of the conduction band (n > 1) the derivative 
changes sign and the mean field increases with increase 
of n. For the number of electrons in the conduction band 



FIG. 6: T^^^^ for the (Ga,Mn)As with Mn concentration 
X = 0.0625 as a function of the electron number n. n = 
corresponds to the system Gao.9375Mno.0625 As with no addi- 
tional donor or acceptor defects. 
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larger than 0.28 per Mn atom the mean field becomes 
positive. Remarkable is the asymmetry of the curve with 
respect to n = 1. The curve is much steeper to the left 
from the point than to the right from it. Correspond- 
ingly, the Curie temperature at n — 2 (one electron in 
the conduction band per one Mn atom) is about 2.5 times 
smaller compared to the Curie temperature at n = (one 
hole in the valence band per one Mn atom) . This reveals 
higher efficency of the holes in the valence band in me- 
diating the exchange interaction between the Mn atoms 
compared with electrons in the conduction band. The 
physical reason for this property is stronger exchange in- 
teraction between the Mn moments and the states at 
the top of the valence band (p-d exchange) compared to 
the exchange interaction between the Mn moments and 
the states at the bottom of the conduction band (s-d ex- 
change) . 

The region of negative n (Fig. ^ corresponds to the 
concentration of of the holes higher than the concentra- 
tion of the Mn atoms. These states of the system can be 
obtained by avoiding the formation of Asca antisites and 
by co-doping with atoms acting as donors. (Of course 
this simple theoretical consideration does not take into 
account the difficulties in the producing of the materials 
with such properties.) 

Moving in the direction of negative n, the value of the 
Curie temperature at first slightly increases. Then the 
mean field value decreases fastly revealing again a trend 
to the antiferromagnetic exchange coupling. The strong 
dependence of the mean field on the number of carriers 
refiects the competition between the antiferromagnetic 
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and ferromagnetic interactions. 

VI. CONCLUSIONS 

We use supercell and frozen-magnon approaches to 
study the dependence of the magnetic interactions in 
(Ga,Mn)As on the Mn concentration. We report the pa- 
rameters of the exchange interaction between Mn spins 
and the MFA and RPA estimates of the Curie temper- 
ature. In agreement with experiment we obtain a non- 
monotonous dependence of the Curie temperature on the 
Mn concentration. We estimate the dependence of the 
Curie temperature on the concentration of the carries in 
the system and show that the decrease of the number of 



holes in the valence band leads to fast decrease of the 
Curie temperature. The strong dependence of the Curie 
temperature on the carrier concentration provides an ex- 
planation to the overestimation of the value of the Curie 
temperature compared to the experiment. 
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